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ABSTRACT 

The Einstein radius plays a central role in lens studies as it characterises the strength 
of gravitational lcnsing. In particular, the distribution of Einstein radii near the upper 
cutoff should probe the probability distribution of the largest mass concentrations in 
the universe. Adopting a triaxial halo model, we compute expected distributions of 
large Einstein radii. To assess the cosmic variance, we generate a number of Monte- 
Carlo realisations of all-sky catalogues of massive clusters. We find that the expected 
largest Einstein radius in the universe is sensitive to parameters characterising the 
cosmological model, especially a& : for a source redshift of unity, they are 

and 54l^ 2 arcseconds (errors denote la cosmic variance), assuming best- fit cosmolog- 
ical parameters of the Wilkinson Microwave Anisotropy Probe five-year (WMAP5), 
three-year (WMAP3) and one-year (WMAP1) data, respectively. These values are 
broadly consistent with current observations given their incompleteness. The mass of 
the largest lens cluster can be as small as ~ 10 15 M Q . For the same source redshift, we 
expect in all-sky - 35 (WMAP5), - 15 (WMAP3), and ~ 150 (WMAP1) clusters that 
have Einstein radii larger than 20". For a larger source redshift of 7, the largest Ein- 
stein radii grow approximately twice as large. Whilst the values of the largest Einstein 
radii are almost unaffected by the level of the primordial non-Gaussianity currently of 
interest, the measurement of the abundance of moderately large lens clusters should 
probe non-Gaussianity competitively with cosmic microwave background experiments, 
but only if other cosmological parameters are well-measured. These semi-analytic pre- 
dictions are based on a rather simple representation of clusters, and hence calibrating 
them with A-body simulations will help to improve the accuracy. We also find that 
these "superlens" clusters constitute a highly biased population. For instance, a sub- 
stantial fraction of these superlens clusters have major axes preferentially aligned with 
the line-of-sight. As a consequence, the projected mass distributions of the clusters 
are rounder by an ellipticity of ~ 0.2 and have ~ 40% — 60% larger concentrations 
compared with typical clusters with similar redshifts and masses. We argue that the 
large concentration measured in A1689 is consistent with our model prediction at the 
1,2a level. A combined analysis of several clusters will be needed to see whether or 
not the observed concentrations conflict with predictions of the flat A-dominated cold 
dark matter model. 

Key words: cosmology: theory — dark matter — galaxies: clusters: general - 
gravitational lensing 



1 INTRODUCTION 

In the standard Cold Dark Matter (CDM) model, which 
for these purposes we shall assume includes the pres- 
ence of a cosmological constant and a flat spatial ge- 
ometry, structure grows hierarchically from small objects 
that merge together to form larger objects (hereafter 
FACDM). Strong gravitational lensing by massive clus- 
ters of galaxies is one of the most important tests of 
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this model in the sense that it probes the rarest high 
density peaks in the universe. For instance, the CDM 
model predicts wide-angle lensing events, on scales as large 
as several tens ar c seconds, due to ma s sive clusters (e.g., 
Turner et alHl984 iNarayan et alJll98i iNaravan fc White! 



19881 : IWambsganss et all Il995l ). This has been broadly 



verified by th e discovery of man y lensed background 
galaxies (e.g., iLe Fevre et al.l Il994l; ILuppino et alJ [19991; 
Gladders et al.ll2003l; IZaritskv fc Gonzalez 2003: Sand et al.l 



2005 



2006) 



Hennawi et al.l [20081 or quasars (jlnada et all 120031 . 
Quantitative comparisons of expected lensing rates in 
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the FACDM model and observed numbers of lenses should 
serve as an important test of our understanding of the uni- 
verse. 

A possible simple test of the FACDM model is the 
statistics of Einstein radii, particularly near the upper cutoff. 
The Einstein radius is a characteristic scale of strong lens- 
ing and is related mainly to the aperture mass it encloses. 
Therefore it is expected that the largest Einstein radii in 
the universe probe the structure and abundance of the most 
massive clusters. This enables a test of the FACDM model 
at the very tail of the halo distribution. An advantage of 
this test is the simple and straightforward determination of 
the Einstein radius in observations and its correspondence 
to identify large lenses. 

Lensing properties of massive clusters have mainly 
been studied using ray-tracing in TV-body simulations 
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1998; iMeneghetti et all 2003a 


. l2005l; iDalal et al. 


2004 


; iHo & White) 120051; Li et al.1 2005 


;lHoresh et al.ll2005l 


Hcnnawi et al. 2007al bl iHilbert et al.l 


2007. 20081). Whilst 



the numerical approach allows one to take account of the 
full complexity of lens potentials, it is often computation- 
ally expensive to perform large box-size simulations which 
retain enough particles in each halo for strong lensing stud- 
ies. In particular reliable predictions for the rarest lensing 
events in the universe require many realisations of such high- 
resolution Hubble-size simulations in order to estimate the 
effect of the cosmic variance. This is impractical with current 
computational capabilities. 

The complementary semi-analytic approaches of- 
ten invoked simple spherically sym metric mass pro- 



files, for calculational re a sons (e.g., Maoz et al. _ 19971: 



Hamana fe Futamasel 1 19971; [Molikawa et al.l 1 19991; ICoorayl 
19991; IWvithe et all l200ll; iTakahashi fc ChiQ | 200l] ; 



Kochanck & White 2001; Kccton 2001; Li & Ostrikcr 2003; 
Lopes fc Milled [20041; iHuterer fc Mai |2004| ; iKuhlen et ail 
20041 ; IChenll2004l 120051 ; lOgurilbOOeT TA more advanced cal- 



culation ad opted an ellipsoid f or proj e cted cluster mass dis- 
tribu t ions (IMeneghetti et al.1 2003e ]; iFedeli fc Bartelmannl 
120071 ; IFedeli et al.1 l2007l. I2008I ). Howev er, because of the 
triaxial nature of FACDM haloes (e. g ., Ijing fc Sutdl2002l ; 
lAlleood et aHl2006l : IShaw et al.ll2006l : lHavashi et al.ll2007h 
the lensing properties of individual clust ers vary drastically 
as a function of viewing angle (e.g., IDalal et al] 120041 ; 



iHennawiet aHl2007hh . resulting in the significant increase 
of average lensing efficiencies due to halo triaxiality. This 
indicates that any analytic models of cluster lensing should 
take proper account of triaxiality for r eliable theoretical 
predictions, as is done in several papers llOguri et al.ll2003l ; 
iQguri fc Keetonll2004l ; [Minor fc Kaplinghatll2008f i~ 

In this paper, we take a semi-analytic approach to pre- 
dict the largest Einstein radius in all-sky survey, based on 
the FACDM model. We invoke an analytic mass function 
of dark haloes to generate a catalogue of massive clusters 
with the Monte-Carlo method. The shape of each halo is 
assumed to be triaxial, and the projection along random di- 
rections is considered. This Monte-Carlo approach allows us 
to evaluate the range of the largest Einstein radii due to 
cosmic variance. We also characterise such "superlens" clus- 
ters, i.e., clusters which produce widest-angle lensing, to see 
how "unusual" are these clusters. 

These issues are well illustrated by detailed observa- 



tions of the largest Einstein radius known to data, which 
may conflict with the FACDM model. The lensing data 
of A1689, one of the best-studied clusters to date, sug- 
gest that t he mass profile is apparently more centrally con- 
centr ated (Broadhu rst et aTll2005al ; iBroadhurst fc Barkanal 
20081) than the FACDM prediction (e. g., iNeto et al.ll2007l : 



iDuffv et al.1 |200S| ; iMaccio' et al] 120081 ), although the ex 
act degree of concentration is somewhat controversial (e.g 
Halkola et alj|2006t iMedezinski et~aill2007l; II imousin et al 



20071 ; IComerford fc NataraiarJl2007r lUmetsu fc Broadhurst 
20081 ). It has been argued that a part of discrepan c y can be 
expla i ned by halo tr i axiality (lOguri et al] l2005al; iGavazzil 
20051 ; iHennawi et ail l2007bl ; ICorless fc Kind 120071 , l2008h 
or the proje ction of the secondary mass peak along the 
line-of-sight (lAndersson fc Madeiskl2004l ; lLokas et al]|2006l : 



iKing fc Corlessl 2007 ), suggesting the importance of careful 
statistical studies with the selection effect taken into ac- 
count. Indeed, it should be pointed out that a weak lens- 
ing analysis of stacked clusters of les ser mass does not ex- 
hibit the high concentrat ion problem l|johnston et al]|2008l ; 
iMandelbaum. et al.ll200sl ). 

We believe that our predictions will be helpful for in- 
terpreting surveys of dista nt (z > 6) galaxies near critical 



curves of massive c lusters (lEllis et alj feOQl; Hu et al.ll2002l; 
Kneib et all 120041 ; iRichard et al] I200II 120081 ; IStark et al.1 



20071 ; 



Willis et al 



120081 : iBouwens et al]|2008h . The survey 



area of this technique is simply proportional to the square 
of the Einstein radius, thus clusters with very large Ein- 
stein radii are thought to be the best sites to conduct this 
search. Our calculations should provide a useful guidance to 
discover such giant lens clusters. 

This paper is organised as follows. We describe our the- 
oretical model in Section [2] Predictions for the largest Ein- 
stein radius and the abundance of large lens clusters are 
shown in Section [3] and Section [4] respectively. Section [5] 
includes the effect of primordial non-Gaussianities. We dis- 
cuss the results in Section [6] and give our conclusion in 
Section [7] Throughout the paper, we consider three cos- 
mological parameter sets obtained from the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP), mainly to show how 
sensitive our results are to cosmological parameters. These 
are the best-fit parameter sets from the WMAP one-year 
data (WMAP1; [Spergel et all 120031 '). (tt M , Q b , fl A , h, n„, 
<t 8 )=(0.270, 0.04 6, 0.730, 0.72, 0.99, .9), WMAP three-year 
data (WMAP3: ISpergel et al.ll2007l ). (0.238, 0.042, 0.762, 
0.732, 0.958, - 761) and WMAP five-year data (WMAP5; 
iDunklev et ail |200S| ), (0.258, 0.044, 0.742, 0.719, 0.963, 
0.796). The most important difference between these models 
is the matter density and the normalisation of matter fluc- 
tuations. Indeed, it has been shown that the smaller values 
of Qm and as in WMAP3 resulted in much smaller num- 
ber of cluster-scale lenses compared with WMAP1 (e.g., 
iLi et al.ll2006ll2007l) . Unless otherwise specified, we adopt 
the WMAP5 cosmology as our fiducial cosmological model. 



2 MONTE-CARLO APPROACH TO THE 
DISTRIBUTION OF EINSTEIN RADII 

We compute the cosmological distribution of Einstein radii 
semi-analytically using a Monte-Carlo technique. First, we 
randomly generate a catalogue of massive dark haloes ac- 
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cording to a mass func t ion. W e assume a fitting formula de- 
rived by I Warren et all (|2006h for the mass function of dark 
haloes: 

0.2538) e^^^^g^) 
; M 2 dlnM v ; 

where M is a halo mass and p(z) is a mean comoving 
matter density at redshift z. We calculate the linear den- 
sity fluctuation o m from the approxima ted transfer func- 
tion presented by lEisenstein fc Hul l|l998h . Throughout the 
paper we adopt the virial mass M — M v i r which is defined 
such that the average density inside a spherical region with 
mass M v i r becomes A(z) times the mean matter density 
of the universe; here A(z) is computed using the spheri- 
cal collapse model in the FACDM universe. For our fidu- 
cial cosmological mod el, A(0) sa 370, A(0.3) « 270, and 
A(l) » 200 (see, e.g., iNakamura fc Sutdfl99fl ). In this pa- 
per we are interested in massive clusters with the masses 
M ~ 1O 15 M0. For comp arison, the ma ss of the Coma clus- 
ter is ~ 1.3 x 10 15 Mm (lHugheslll989l), and tha t of A1689 
is ~ 2.1 x 1O 15 M l|llmetsu fc Broadhurstll2008l '). With this 
mass function, the number of dark haloes for each redshift 
and mass bin can be written as 

N = 3-777 AzAM = SlD A {zf — ( 
dzdM az 

with Da(z) and dr/dz being the angular diameter distance 
and the proper differential distance, respectively. Through- 
out the paper we adopt the solid angle of Q — 40, 000 deg 2 
which roughly corresponds to all-sky excluding the Galac- 
tic plane. A realisation of dark haloes is then constructed 
by computing the expected mean number of dark haloes for 
each bin adopting Az = 0.01 and A(log M) — 0.02, and 
generating an integer number from the Poisson distribution 
with the mean. The Monte-Carlo catalogues are generated 
in the range of cluster masses larger than the minimum 
mass M min . We adopt M mln = 4 x 10 14 M Q for WMAP1 
and Mnun = 2 x 10 14 Af Q for WMAP3 and WMAP5, which 
are sufficiently small not to affect our results. On the other 
hand, the maximum cluster masses for these cosmologies are 
A/max ~ 3 - 5 x 1O 15 M (see 

Each dark halo is assum ed to have a triaxial shape. 
Following Ijing fc Sutcl (|2002l . hereafter JS02), we model the 
density profile as 



r [Mpc] 
0.1 1 



10 



zf^AzAM, (2) 



p(R) 



R 



Scdpcrit(z) 



1 



(R/Ro)(l + R/Ro) 2 I 1 + (R/Rt) 



x 2 v 2 z 2 



(3) 



(4) 



The m odel is a triaxial generalisation of the iNavarro et al.l 
1 19971 . hereafter NFW) density profile. The concentration 
parameter for this triaxial model is defined by c e = 
Re/Ro, where R e is determined such that the mean den- 
sity within the ellipsoid of the major axis radius R e is 
5A(z) (c 2 /afc) p(z)(l+z) 3 . The characteristic density <5 c d 
is then written in terms of the concentration parameter. 
As suggested in JS02 we relate R e to the virial mass M vlr 
of the halo by adopting a relation i? e /r v i r = 0.45, where 
r v i r is spherical virial radius computed from the virial mass. 
A change from JS02 is the inclusion of a truncation term, 
[1 + (R/R t ) 2 ]~ 2 , such that the r adial profile does not ex- 
tend far beyond the virial radius (|Baltz et al.ll2008l see also 
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Figure 1. The illustration of various characteristic scales for 
large clust er lenses. The best-fit spher ical NFW density profile 
of A 1689 llUmetsu fc Broadhurstl |2008|) p(r) normalised by the 
mean matter density of the universe at that redshift p is plot- 
ted as a function of the radius. Thick and thin lines show the 
NFW profiles with and without the truncation term (see eq. [3]). 
The Einstein radius for the source redshift of unity, the scale ra- 
dius of the spherical NFW profile, and the virial radius of the 
cluster are indicated by vertical dotted lines with labels of rg, 
r s (~ J?o), and r v j r , respectively. The total mass of a sphere de- 
fined by the Einstein radius is 1.1 X 10 14 AfQ. and the cylindrical 
mass projected within the Einstein radius is 2.0 X 1O 14 M0, which 
should be compared with the virial mass (total mass inside r v ; r ), 
M vir = 2.1 X 10 15 M Q . 



iTakada fc Jainl |2003h ■ We choose R t = 4r v i r which can be 
translated into the truncation at roughly twice the virial 
radius for massive haloes. We note that the truncation is 
introduced for haloes with very small c e ; in this situation 
mas ses outside the virial r adii dominate the lens potentials 
(see lOguri fc Keetonll2004t) , and thus the truncation is nec- 
essary to avoid such unrealistic situations. The truncation 
has a negligible effect on the Einstein radii of most haloes. 

To give a rough idea of various length scales for mas- 
sive lensing clusters, in Figure [1] we plot the best-fit ra- 
dial NFW density profile o f A1689 derived from lensing 
|Umetsu fc Bro adhurst 2008). The Einstein radii of massive 
lensing clusters are typically ~ 5% of the virial radii r v j r 
which are a few Mpc for these clusters. The density at the 
Einstein radius is ~ 10 5 times more than the mean matter 
density of the universe p. The Figure also indicates that our 
truncation of the NFW profile (see eq. [3]) only affects the 
radial density profile at r > r v i r - The radial profile crosses 
the mean matter density p at several times the virial radius. 

The axis ratio and concentration parameter for each 
halo are randomly assigned according to the probability dis- 
tribution functions (PDFs) derived by JS02. Specifically the 
probability distributions for triaxial axis ratios are given by 
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p(a/c) 



p(a/b\a/c) 



27T(T, 



exp 



(r ac - 0.54)^ 



2a] 



dr* 



d(a/c) ' 



2(1 - r min 



0.07[f2(z)]° 



2a /b — 1 — r m i r , 

1 1 min 



_ a ( M \ 
= Max[o/c, 0.5] 



(5) 
(6) 

(7) 
(8) 



The characteristic nonlinear mass M* is defined such that 
the overdensity at that mass scale becomes 8 C = 1.68. The 
best-fit value for the width of the axis ratio distribution a s 
is o s = 0.113. Note that p(a/b\a/c) — for a/b < r m i n . On 
the other hand, the probability distribution for the concen- 
tration parameter is well approximated by the log-normal 
distribution: 



2ll(T, 



exp 



(In c e — In c e 
2aJ 



(9) 



with the width of the distribution a c — 0.3. We include a cor- 
relation between the axis ratio and concentration parameter 
by adopting the fo llowing form for th e median concentration 
parameter c e (see lOguri et aljfeoOSl ): 



Max[/ C ,0.3]A E 



A(z c 



f c — 1.35 exp 



A( 
03^ 2 



1 + Zc 



3/2 



(10) 



(11) 



where A e = 1.1 for FACDM. The collapse redshift z c of 
a halo with mass M is estimated by solving the following 
equation involving the complementary error function: 



erfc- 



S c (z) - Sg(0) 



fM 



(12) 



Here the linear overdensity at redshift z is S c (z) = S c /D(z) 
with D(z) being the linear growth rate, and ct/m an d &m 
are linear density fluctuations for the mass scales of fM 
and M, respectively. Note that a 2 is computed at z — in 
this equation. We adopt / = 0.01 following JS02. Equation 
(|11|) suggests that the concentration parameter becomes too 
small for very elongated haloes (a/c <C 1). Since the fitting 
formula of f c was derived at f c > 0.3 (see JS02), we modified 
the prefactor in equation (|10|) from f c to Max[/ C , 0.3] in or- 
der to avoid unrealistically small values of the concentration 
parameter. 

We need to specify the orientation of each halo relative 
to the line-of-sight direction to compute lensing properties. 
The orientation of dark haloes can be specified by the follow- 
ing two angles; a (0 < a < n) defined by an angle between 
the major axis of the triaxial halo and the line-of-sight direc- 
tion, and [3 (0 < j3 < 2n) which represents a ro tation angle 
in a plane perpendicular to the major axis (see lOguri et al.l 
l2003l ; lo"guri fc K ccton 20Q3)- Assuming that the orientation 
of each halo is random, the PDFs of these angles are given 
by 



p(a) = 



(13) 
(14) 



We perform the pr ojection of the tria xial halo following the 
procedure given bv lQguri et ail (|2003T ) and compute the pro- 
jected convergence and shear maps for a given source red- 
shift z s : 



. . &TNFW j, 

H*,!/) = — 7, — /nfw 



1 x 2 y 2 



&TNFW = 



45 cd p crit (z)Ro 



VfXc 



dz- 



{l + (x 2 + z 2 )/x 2 Y 



(15) 



(16) 



(17) 



/NFw(#) — . - ______ _______ 

' Vx 2 + z 2 (l + Vx 2 + z 2 ) 2 ' 

where E cr is the critical surface mass density for lensing and 
xt — 4r v i r /i?o is the truncation radius. The parameters q x , 
q y , and / are complicated functions of the ax is ratios and the 
projection direction (see lOguri et al.l 120031 . for explicit ex- 
pressions). The axis ratio of the projected mass distribution 
is given by q = q y /q x - This model introduces the ellipticity 
in the projected density, and therefore does not suffer from 
unphysical mass distributions that are see n if the elliptic- 
ity is introduced in the lens potential (e.g.. lGolse fc Kneibl 
2002). Since critical curves of projected triaxial haloes are in 
general neither circles nor ellipses, the definition of the Ein- 
stein radii for these systems are not trivial. In this paper we 
compute the Einstein radii as follows. First we calculate dis- 
tances from the halo centre to the (outer) critical curve along 
the major and minor axes of projected two-dimensional den- 
sity distribution, which we denote 9 X and 9 y , respectively. 
Then we estimate the Einstein radius of the system by the 
geometric mean of these two distances: 



(18) 



By computing Einstein radii for all the massive dark haloes 
we have randomly generated, we obtain a mock all-sky cat- 
alogue of Einstein radii. For each model we consider, we 
generate 300 of all-sky realisations in order to assess the 
cosmic variance of the largest Einstein radii in the universe. 

To demonstrate the importance of triaxiality on this 
study, we compute Einstein radii of 1,000 massive dark 
haloes with the virial mass M vir = 2 x 1O 15 M and the 
redshift zi — 0.3. The concentration parameter, axis ratios, 
and the orientation with respect to the line-of-sight direc- 
tion of each halo are randomly generated using the PDFs 
described above. We show the resulting distribution of 8e in 
Figure [5] The Figure indicates that haloes of the same mass 
can have a wide range of the Einstein radii. They are corre- 
lated with the orientation of the halo such that largest Ein- 
stein radii are caused only when the major axis of haloes is 
almost aligned with the line-of-sight direction (| cosa| ~ 1), 
implying a strong orientation bias in large lens clusters. 



3 LARGEST EINSTEIN RADIUS AND 

PROPERTIES OF THE LENSING CLUSTER 

3.1 Probability distribution of the largest 
Einstein radius 

First we take the cluster that has the largest Einstein ra- 
dius from each all-sky realisation. From the 300 realisations 
for each cosmological model, we can not only construct a 
probability distribution of the largest Einstein radius in the 
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Figure 2. The distribution of the Einstein radii $e for 1,000 



triaxial haloes with fixed virial mass M v 



2 X lO 15 A/ . The 



haloes are located at zi = 0.3, and the source redshift is assumed 
to be z s = 1. The distribution is shown as a function of | cosa|, 
where a is an angle between the major axis of each halo and the 
line-of-sight direction. For comparison, the virial radius of the 
cluster corresponds to 9 ~ 640". 



universe, but also obtain expected properties of the lensing 
cluster. In what follows, we consider three source redshifts, 
z 3 = 1, 3, and 7. The source redshift of z s = 1 is more rele- 
vant to typical giant luminous arcs or weak lensing studies, 
whereas results for z 3 — 7 are more important in searching 
for high-redshift galaxies near critical curves. The redshift 
distribution of strongly-lensed faint backgr ound galaxies in 
massiv e clusters has a peak at z s ~ 3 fe.g.. lBroadhurst et al.l 
2005b) . The probability distributions of the largest Einstein 
radius in all-sky, # max , are shown in Figure [3] As expected, 
#max is quite dependent on the cosmological model. For 
the source redshift z s = 1, the median of the largest Ein- 
stein radius is max = 54" for WMAP1, max = 35" for 
WMAP3, and max = 42" for WMAP5. The different values 
of as (a 8 = 0.9 for WMAP1, erg = 0.76 for WMAP3, and 
erg = 0.8 for WMAP5) change the abundance of massive 
dark haloes and its redshift evolution drastically, which is 
why the largest Einstein radius is sensitive to <7g. It is also 
found that the value of the largest Einstein radius is quite 
dependent on the source redshift as well: # max for z a = 7 is 
approximately twice as large as that for z s — 1. 

The cluster which has the largest known Einstein 
radius to date is A1689. It is a massive cluster lo- 
cated at low redshift (zi = 0.18), and its Einstein ra- 
dius is well constrained from many multiply-lensed back- 
ground galax ies and weak lensing to be 8e = 45" 
for z a = 1 dTvson et al I |l99fJ: iMiralda-Escude fc Babull 



19951; IClowe fc Schneiderl l200l| ; [Broadhurst et al.l l2005al lbf 
Broadhurst 2008). From the best- fit mass model 
of Umetsu fc Broadhurst! l|2008l ). we derive #e = 53" for 
z a = 3 and #e = 56" for z s = 7. We emphasise that these 
values should be viewed as the lower limits of # max , since 
even larger lens clusters may be discovered in the future. 
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Figure 3. Probability distributions of the largest Einstein radius 
9 maI , constructed from 300 Monte-Carlo realisations of all-sky 
massive cluster catalogues. From top to bottom panels, source 
redshifts are assumed to be z s = 1, 3, and 7. Results for three 
different cosmological models, WMAP1 (dashed), WMAP3 (dot- 
ted), and WMAP5 (solid) are shown. Arrows indicate the values 
of the Einstein radii of A1689, which have the largest known Ein- 
stein radii. Note that they correspond to the lower limits of # max 
in observations (see text for details). 



Nevertheless, we compare these values with our theoretical 
predictions in Figure [3] For z s — 1, the Einstein radius of 
A1689 is quite consistent with our prediction for our fidu- 
cial cosmological model, WMAP5. On the other hand it is 
slightly larger (smaller) than our prediction for WMAP3 
(WMAP1), but is within 90 percentile for both WMAP1 
and WMAP3. In contrast, the observed values are consis- 
tent with WMAP3 and lower than WMAP1 and WMAP5 
for Zs — 3. For the higher source redshift z s = 7, the Ein- 
stein radius of A1689 is even smaller than our prediction for 
WMAP3. An implication of this is that there may exist a 
cluster with the Einstein radius larger than that of A1689 
for that source redshift. We note that cosmological models 
which predict 9 mjI close to the Einstein radius of A1689 for 
z 3 = 7, such as models with as even smaller than WMAP3, 
should predict too small # max for z s = 1 to be consistent 
with A1689. 

Whilst the better match to WMAP3/WMAP5 is con- 



© RAS, MNRAS OOClTHTTl 



6 M. Oguri and R. D. Blandford 



1.5 - 



0.5 - 



71 1 1 — I I I. I I 

WMAP5/ " 




10 15 
M vir [M ] 



10 16 



Figure 4. The distribution of the mass and redshift of the clus- 
ter producing the largest Einstein radius. Each point corresponds 
to one Monte-Carlo realisation. Results are shown for three dif- 
ferent redshifts, z s = 1 (pluses), 3 (open squares), 7 (crosses). 
The WMAP5 cosmology is assumed in this plot. Also plotted 
are contours of constant X-ray fluxes, inferred from the correla- 
tion between bolome tric X-ray luminosities and halo virial masses 
llShimizu et al 1 l2003l, assuming no redshift evolution). From right 
to left, the contours indicate X-ray fluxes of fx = 10 -11 , 1CT 12 , 
10" 13 , and 10" 14 ergs^cm- 2 . 



sistent with recent results from th e cluster abundance (e.g., 
iDahlel l200fj ; iGladders et al] 120071 ; iMantz et all l2008h , the 
large discrepancy between the predicted Einstein radii for 
WMAP1 and those of A1689 does not necessarily exclude 
WMAP1 cosmology as the current observed f? max correspond 
to the lower limits. Complete surveys of large lens clus- 
ters are necessary to extract useful cosmological information 
from this statistics (see i|6.ip . 



3.2 Cluster mass and redshift 

Next we examine the expected mass and redshift distri- 
bution of the cluster which produces the largest Einstein 
radii. The distributions shown in Figure|4]indicate that wide 
ranges of mass and redshift are possible. In particular, it is 
worth noting that the mass of the cluster can be as small 
as M v i r < 1O 15 M0. In addition, the cluster can be located 
at quite high-redshifts, up to zi ~ 1 and beyond, in the 
case of z 3 — 3 and 7, which will be difficult to access via 
X-ray observations with currently operating telescopes (see 
also discussions in ij6.1[) . On the other hand, for z s = 1 the 
lens cluster is likely to be located at zi < 0.5. Clearly the 
diversity of the mass and redshift is a consequence of halo 
triaxiality, which we will explore later. 

Here a natural question to ask is whether or not the 
cluster having the largest Einstein radius is the most mas- 
sive cluster in the universe. To check this we take the most 
massive cluster in each realisation and construct the proba- 
bility distribution of its mass. It is then compared with the 
PDF of the mass of the largest lens cluster. We show the re- 
sult in Figure [S] As we discussed, the cluster with # max has 
a wider range of the mass and thus is not necessarily the 
most massive cluster in the universe. However, the overlap 
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Figure 5. Probability distributions of masses M v i r of the most 



massive cluster in the universe (thick solid) and the cluster pro- 
ducing the largest Einstein radius (thin solid for z a = 1 and 
dashed for z s = 7). From top to bottom panels, we show results for 
the WMAP1, WMAP3, and WMAP5 cosmologies, respectively. 



of the PDFs at the high-mass end for all the three cosmo- 
logical models suggests that in some cases the largest lens 
corresponds to the most massive cluster. It is interesting to 
note that the PDFs for z s — 7 extend to lower masses than 
those for z s = 1. This is because the most massive clusters 
are typically located at zi ~ 0.1 — 0.4, whereas the geomet- 
rical lensing efficiency for the source redshift z s = 7 is the 
highest at around zi ~ 1 where clusters are on average less 
massive (see also Figure Q. 

3.3 Expected properties of the lensing cluster 

It has been argued that the population of lenses is markedly 
different from that of nonlcnscs in several ways (e.g. , 



Oguri et al.|2005bl;lHennawi et alj|2 007b; M oller et ai1l2007l ; 



Fedeli et all 120071 ; iRozo et al.ll2008bl ). The largest Einstein 
radius represents the most extreme case of lensing clusters, 
which suggests that the cluster population may be biased 
even more strongly. Here we quantify the lensing bias in the 
cluster producing r? max from our Monte-Carlo realisations. 

As discussed above, an important parameter here is 
the orientation of the cluster, specifically the angle a be- 
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Figure 6. Probability distributions of the triaxial concentration parameter c e (top), the triaxial axis ratio a/c (centre), and the angle 
between the major axis and the line-of-sight direction a (bottom) for the cluster producing the largest Einstein radius are plotted. For 
comparison, the PDFs for typical (unbiased) clusters with the same mass and redshift distributions are shown by dashed histograms. 
From left to right panels, we change the source redshift from 1 to 7. Although results shown here are for the WMAP5 cosmology, we 
confirmed that the PDFs for the other two cosmologies arc similar. 



tween the major axis and the line-of-sight direction. An- 
other important parameter is the minor-to-major axis ratio, 
a/c, since the projection effect is stronger for more triaxial 
clusters. Finally the concentration parameter of the triax- 
ial model, c e , should also be of interest because the strong 
lensing efficiency is known to be sensitive to the halo con- 
centration. 

Figure [U] shows probability distributions of these three 
parameters for the cluster having # max . We also show the 
PDFs for a "typical" cluster which has the same mass and 
redshift probability distributions as those of the lens cluster 
but the axis ratio, the concentration, and the orientation are 
re-assigned from their original PDFs. Thus the comparison 



of the lens and typical cluster PDFs provides the degree of 
lensing biases. Strikingly, we find that the cluster is almost 
always aligned with the line-of-sight direction. For instance, 
the probability of having | cos a\ > 0.9 (a < 25.8°) is 0.88 for 
z s = 1 and 0.95 for z s — 7 for the WMAP5 cosmology. It is 
also found that the lens cluster is more triaxial than typical 
clusters. Because of the correlation between c e and a/c, the 
triaxial concentration parameter for the lens cluster becomes 
smaller. The strong biases in the orientation and triaxiality 
indicate that such largest lens cluster is indeed very unusual 
in terms of its internal structure and configuration. 

The projection effect of the triaxial halo has a large 
impact on the apparent mass profile constrained from lens- 
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Figure 7. Similar to Figure[6] but the PDFs of the projected two-dimensional virial mass M2D, the projected two-dimensional spherical 
concentration parameter C2D, and the projected ellipticity e2D are shown. 



ing data (| Oguri et al, 2005a). To investigate this, we char- 
acterise the projected two-dimensional (2D) mass distribu- 
tion of each triaxial cluster by the following three param- 
eters: the mass M2D and concentration parameter C2D of 
the NFW profile, which are obtained from the projected 
surface mass density by ignoring the elongation along the 
line-of-sight, and the ellipticity of the surface mass density 
62D- Our parameter C2D corresponds to the standard con- 
centration parameter which has been studied from analysis 
of observed lensing clusters, and thus is useful for discussing 
possi ble high concentrations f rom lens mass reconstructions 
fe.g.. lBroadhurst et al-lkoOSal ). We can relate these param- 
eters to those of the triaxial halo by simply comparing the 
expressions of the projected surface mass densities (K. Taka- 
hashi et al., in preparation): 



r(M 2 ; 



= RoQx 



&NF\v(Af2D, C2D 



&TNFw(A/vi r , C e 



(20) 
(21) 



e2D = q, 



(19) 



where 6nfw is the lensing strength parameter for the spher- 
ical NFW profile. 

We show probability distributions of these 2D parame- 
ters, for both the cluster having # ma x and corresponding typ- 
ical cluster, in Figure[7] It is clear that the largest lens cluster 
is highly biased in terms of the 2D parameters as well. In 
lensing observations the lens cluster looks more massive and 
more centrally concentrated than typical clusters. Indeed 
this is expected from the strong orientation bias (see above), 
because both the mass and concentration of the cluster pro- 
jected along t he major axis are known to be significa ntly 
overestimated (jOguri et alj|2005al : ICorless fc Kingj|200St ). In 
addition we find that the cluster should appear rounder. 
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Table 1. The largest Einstein radii and expected properties of the lensing clusters. We show median values and 68% confidence intervals 
estimated from 300 Monte-Carlo realisations. Values in parentheses indicate corresponding parameter values for typical clusters: they 
are estimated by adopting same mass and redshift distributions as those of the lens cluster and re-assigning concentrations, axis ratios, 
and orientations according to the PDFs. 
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Its expected median ellipticity of 0.25 — 0.3 is significantly 
smaller than that of typical triaxial cluster, 0.45—0.5. Again, 
this is because of the orientation bias: the projected mass 
distribution of a triaxial halo is most elliptical when it is 
projected along the middle axis, and least elliptical when 
projected along major or minor axis. This means that the 
circularity of projected density distribution do not necessar- 
ily imply the sphericity of the cluster. We summarise our 
numerical results in Table [T] 

We are now in a position of discussing the high 
concentration parameters observed in some lens-rich clus- 
ters. For instance, from the strong and weak lensing data 
lUmetsu fc BroadhurstJ l|2008l ) refined the concentration pa- 
rameter of A1689, which has the largest known Einstein ra- 
dius, to be C2D = 12.7 assuming a spherical NFW profileQ 
Adopting the source redshift z s = 1, our model predicts 
that such a cluster should have the concentration parame- 
ter in the range 6.7 < C2D < 11.8 at 68% confidence and 
5.5 < C2D < 14.8 at 90% confidence. Thus we conclude that 
the high concentration of A1689 is consistent with the the- 
oretical expectation based on the FACDM model at 1.2<r 
level. We note that the redshift of the cluster is also con- 
sistent at ~ la level. However the halo concentration of 
A1689 is still slightly larger than the theoretical expectation. 
Therefore statistical studies of concentrations for several 



1 Although iLimousin et al.l d2007f) obtained somewhat smaller 
value of concentration, C2D = 9.6, form their strong and 
weak lensing analysis, it has been argued that their best-fit 
model pre dicts the Einstein radius sm aller than what is ob- 
served (sec Umctsu & Broadhurst 2008). Since the Einstein ra- 
di us is of central interes t , in this paper we adopt the result 
of lUmetsu fc BroadhurstJ (1200811 as a fiducial best-fit model of 
A1689. 



large lens clusters, as attempted in iBroadhurst fc Barkanal 
( 2008) and IBroadhurst et ail (|2008l ). will be essential to as- 
sess whether or not the large concentration problem is in- 
deed existent. 



4 DISTRIBUTION OF EINSTEIN RADII 

Thus far we focused our attention on the largest Einstein 
radius on the sky. Another interesting quantity to investigate 
is the number of clusters which have relatively large Einstein 
radii. Here we derive the expected number distribution of 
such large Einstein radii from our Monte-Carlo realisations. 

In Figure [S] we plot the cumulative number distribu- 
tions for all the three cosmological models. It is found that 
the number decreases exponentially with increasing Einstein 
radius. As in the case of the largest Einstein radius, the 
abundance of large lens clusters is quite sensitive to the cos- 
mological model. For instance, the all-sky numbers of clus- 
ters which have 8e > 20" (for z s = 1) are predicted to be 
~ 150, ~ 15, ~ 35, for WMAP1, WMAP3, and WMAP5 cos- 
mologies, respectively. The large difference of the cumulative 
numbers between WMAP1 an d WMAP3 is broadly consis- 
tent with lLi et ail (|2006lJ2007h who investigated strong lens- 
ing probabilities in two different cosmological simulations. 
The result suggests that it provides useful constraints on 
cosmological parameters. 

One of the main findings about the largest Einstein ra- 
dius was that the lens clusters constitute a highly biased 
population (see Section 13.3 P . One might expect that the 
strong bias is due to the rareness of such largest lens, and 
thus it is interesting to check lensing biases for more com- 
mon lens events. In Figure [9] we show biases in clusters 
which have Einstein radii larger than certain values. Here 
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Figure 8. Expected all-sky cumulative number distributions of 
large Einstein radii. The median and lcr error-bars of the numbers 
are estimated from 300 Monte-Carlo realisations. Results are pre- 
sented for all the three cosmological models, WMAP1 (crosses), 
WMAP3 (squares), and WMAP5 (triangles). From top to bottom 
panels, results for z s = 1, 3, and 7 are shown. 



we adopt two parameters to quantify lensing biases. One is a 
measure of the orientation bias, f\ cos Q >o.9 , which is defined 
by the fraction that the angle between the major axis of the 
cluster and the line-of-sight direction satisfies | cosq| > 0.9. 
If the orientation of clusters is completely random, we should 
have /| COS a>o.9| ~ 0.1. The other parameter is to describe 
the 2D concentration bias, C2rJ/c2D,ran, which is the ratio of 
median 2D concentration parameters (see also Section [3,3p 
among large lens clusters and corresponding typical clusters 
with similar mass and redshift distributions. The parameter 
becomes unity if no lensing bias is present. 

The result shown in Figure[9]indicates that clusters with 
large Einstein radii are similarly highly biased as the cluster 
producing the largest Einstein radii in all-sky. We find that 
the degree of the orientation bias depends on the limiting 
Einstein radius. Populations of clusters with smaller limiting 



Einstein radius, which are regarded to represent less extreme 
populations of lensing clusters, are less biased in terms of 
their orientations. On the other hand, the degree of the 2D 
concentration bias does not depend strongly on the limiting 
Einstein radius. The bias in the 2D concentration parame- 
ter comes both from the enhancement of the apparent con- 
centration due to the orientation bias and from the bias in 
the 3D triaxial concentration parameter. The enhancement 
of C2D due to the orientation bias is larger for more triax- 
ial haloes, however because of the anti-correlation be tween 
the axis ratio and concentration (|jing fc Sutoll2002r i such 
haloes are intrinsically less concentrated (i.e., have smaller 
c e ). The behavior is therefore expected to reflect the com- 
plicated combination of these two biases which are more or 
less counteract with each other. 

Our resu l ts can be compared with those of 
iHennawi et~aH (|2007bh who analysed lensing biases using 
ray-tracing of iV-body simulations. Their qualitative results 
are similar to ours: they found that lensing clusters tend 
to have the major axis aligned with the line-of-sight and 
larger 2D concentrations. However, the quantitative results 
are different. Their orientation bias of /| CO sa>o.9| ~ 0.25 
and 2D concentration bias of c^u/c2u,ia,n ~ 1-34 are smaller 
than our results (see Figure |5J). We ascribe this difference 
to the differen t definit ions of the lens cluster populations. 
IHennawi etall (|2007bl ) derived the distributions for lens 
clusters by calculating those from all clusters with a 
weight of lensing cross sections, without any restriction 
to the Einstein radii. Therefore their results are relevant 
to more common lens clusters with smaller Einstein radii, 
say 10" — 15", whereas our results are applicable only to 
superlens clusters with unusually large Einstein radii. Our 
finding that the orientation bias decreases with decreasing 
Einstein radius is consistent with this interpretation. 



5 PRIMORDIAL NON-GAUSSIANITY 

The results presented so far are based on standard uni- 
verses evolved from Gaussian initial conditions. Since the 
abundance of massive clusters and its redshift evolution are 
know n to be very sensitive to primordial non - Gaussianities 



(e.g.. iMatarrese et al 



20001: IVerde et al l 
120071 ; ISadeh et al.l 



2001 



iMathis et~al 



2007: Fedeli et al 



2004: IGrossi et al . 

20081 : iDalal et al.l l2008f T our statistics are also expected 
to be dependent on primordial non-Gaussianities. The ef- 
fect of primordial non-Gaussianities is particularly of im- 
portance given possible detections in th e cosmic microwave 
background (CMB) an isotropies (e.g., IVielva et al. 1 12004 
lYadav fcW andelt 2008j). In this section, we repeat the same 
calculations conducted in the previous sections, but includ- 
ing levels of primordial non-Gaussianities currently of inter- 
est. 

In order to quantify the effect of primordial non- 
Gaussianities, in this paper we ad opt the local non- 
Gaussianity of the follow ing form fe.g.. [Komatsu fc Spergell 
l200ll : iBartolo et al1l2004r i : 



$ = 



0+/NL (0 2 -<0 2 >) 



(22) 



where <E> is the curvature perturbation and (f> is an auxil- 
iary random- Gaussian field. The level of primordial non- 
Gaussianity is characterised by /nl, which we assume con- 



© RAS, MNRAS 000.ITHT71 



The largest Einstein radius 11 






1 




0.8 


OJ 




o 

A 


0.6 


8 




m 
O 
o 


0.4 








0.2 









I I I I I I I I I I I I I I I I I I I 



1.8 

1 1-6 

Q 
CM 

o 1.4 
IcP 1.2 
1 



20 30 40 50 
E [arcsec] 
i i i i i i i i i i i i i i 



40 50 60 70 
E [arcsec] 



40 50 60 70 80 
E [arcsec] 



z =1 



I I I I I I I I I I I I I I I I I M 

1.8 F 4 1.8 

z_=3 



^^*" A *"*^Sn itxxxxx x xxxxx 
xX xxxxxxx x x>A*<xxx xx 



I I I I I I I I I I I I I I 



2 1-6 
d 

CM 

« 1.4 
lcFl.2 



I I I I I I I I I I I I I I I I I I I I 



£ 1-6 

Q 

CM 

o 1.4 



^Ic?1.2 



1 



_l 1 1 1 1 1 


II | II II | II M_ 










— m. - ^ 






1 1 1 1 1 1 


i i 1 i i i i 1 i i i i 



20 30 40 50 
E [arcsec] 



40 50 60 70 
6 E [arcsec] 



40 50 60 70 80 
E [arcsec] 
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the case with no bias. Crosses, open squares, and filled triangles denote results for WMAP1, WMAP3, and WMAP5 models, respectively. 
We consider the source redshift of z s = 1 (left), 3 (centre), and 7 (right). 



stant. In this model positive /nl corresponds to posi- 
tive skewness of the density field. In observation, CMB 
anisotropies h ave constrained i ts val ue to be |/nl| < 0(100). 
For instance, iKomatsu et aL | (|2003l ) derived -58 < /nl < 
134 at 95% confidence from the tempe rature map of the 
WMA P first-year data. More recently, lYadav fc Wandeltl 
(2008) claimed the detection of positive /nl, 26.9 < 
/nl < 146.7 at 95% confidence, from the improved anal- 
ysis of the WMAP t hree-year data. On the other hand, 
lKomats"u" et al.l (|2008h claimed that the WMAP five -year 
data are marginally consistent with Gaussian fluctuations 
at 95% confide nce, —9 < /nl < 1 11. 

We follow iDalal et"al] (|2008l ) to calculate the number 
of massive haloes in this non-Gaussian model. From ana- 
lytic considerations and A^-body simulations, they derived a 
simple fitting formula of the mass function: 



rfWNG 

dM 



dM dn 



M dM J2^a 



■ exp 



7 



(M/M - fuf 
2a) 



.(23) 



where dn/dMo is a halo mass function with a Gaussian ini- 
tial condition, which we adopt equation (UJ, and 

(24) 

a f = 1.4 X 1U -(J/NL|<T 8 r -<T Mo - (25) 

From this expression, we can see that the positive /nl results 
in the enhancement of the abundance of massive clusters. 
Primordial non-Gaussianities affect not only the abun- 



fii = 1 + 1.3 x 10 4 /nl0- 8 o- m ^, 
1.4 x IO^G/nlM"' 8 ^ 1 



dance of massive clusters but also their formation histories. 
Since the concentrations of d ark haloes are correla ted with 
their mass assembly histories (We chsler et al.l2002l ). primor- 
dial non-Gaussianities should have an impact on the halo 
concentration parameter as well . Indee d, A^-body simula- 
tions done by lAvila-Reese et all (|2003l ) showed that posi- 
tive (negative) skewness in the initial density field results 
in larger (smaller) halo concentrations. We crudely include 
this effect by modifying the linear ov erdensity in equation 
O as follows l|Matarrese et al.ll2000l ): 



dc(z) -» S c (z)y/l-S 3 Sc(z)/3. 



(26) 



We estimate t he skewness 53 as S3 ~ 6/nl o M a^, with 
4 x 10~ 5 (|Dalal et al.l l2008h. In this model, primordial 
non-Gaussianity of /nl = ±200 translates into the modifi- 
cation of median concentration parameters for haloes with 
the mass 10 15 Mq by ~ ±5%. Give n the current level of 
constraints on /nl fr om WMAP (e.gjKomatsu et al.ll2003l; 
ISpergel et afl 120071 ; lYadav fc Wandeltl 120081 : iHikage et all 
2008), in this section we consider three non-Gaussian mod- 
els, /nl = —100, 100, and 200, as well as the Gaussian case 
/nl = studied in the previous sections. For each model, we 
compute 300 realisations of all-sky cluster catalogue to esti- 
mate median and cosmic variance of large lenses. First, we 
examine the effect of primordial non-Gaussianities on the 
probability distribution of the largest Einstein radii # max . 
The dependence of # max on /nl is displayed in Figure 1101 
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We find that primordial non-Gaussianity of |/nl| ~ 100 
hardly affect S max . Although # max is increased by ~ 10% 
from /nl = —100 to 200, it is clearly smaller than the cos- 
mic variance. This suggests that the observation of # max 
hardly constrains /nl, at least not so tightly as the cur- 
rent WMAP data do. On the other hand, the plot shown in 
Figure[TT]suggests that we can in principle detect primordial 
non-Gaussianities of |/nl| ~ 100 from the all-sky abundance 
of clusters with relatively large Einstein radii, N(> 9e)- For 
instance, iV(> 20") for z 3 = 1 is 34 ± 6 and 60±g° (68% 
confidence) for /nl = and 200, respectively. The abun- 
dances of large lenses for higher source redshifts are more 
sensitive to /nl, because lens clusters are located at higher 
redshifts where the cluster abundance is more sensitive to 
primordial non-Gaussianities. One of the reasons why the 
abundance of large Einstein radii is better in probing /nl 
than the observation of 9 mal is its smaller cosmic variance. 

The results presented above suggest that constraints 
on /nl are not improved very much by including large 
lenses, compared with current CMB constraints. How- 
ever, some inflation models predict st rongly scale-dependen t 
primordial non-Gaussianities (e.g., lLoVerde et al.l [2008), 
and thus independent constraints from clusters of galax- 
ies, which probe smaller scales (a few Mpc) than CMB 
anisotropics (> 100 Mpc), can be very important to test 
such scale-dependence. The best constraints on primordial 
non-Gaussianities at the cluster scale are expected to be 
obtained by the number count of clusters at high-redshifts, 
detected in radio, X-ray, or optical, but an accurat e calibra- 
tion of cluster masses is alw ays challenging (e.g., iHu et all 
120071 ; iTakada fc Bridle! I2007T ) . The statistics of large lenses 
may therefore provide an important complementary test of 
cluster-scale primordial non-Gaussianities. 



6 DISCUSSIONS 

6.1 Observational strategy 

In this paper, we derived all-sky distributions of large Ein- 
stein radii based on the FACDM model. An advantage of 
the Einstein radius statistics is that it is determined quite 
well from observations, provided that strongly lensed arcs 
are ob served. How many arcs do we expect? lOguri et al.l 
(2003) computed the number of lensed arcs in a typical 
massive cluster, with a mass of 2 - 3 x 1O 15 M , to be 
~ 1 for the arc magnitude limit of ~ 26 mag. Since the 
clusters studied in this paper have 2 — 3 times larger Ein- 
stein radius than typical clusters of similar masses, the 
expected number of lensed arcs for these clusters should 
also be larger by a factor of 5 — 10. Therefore, we con- 
clude that reasonably deep (~ 26 mag) optical imaging 
of clusters having large Einstein radii should always re- 
veal several strongly lensed arcs, which will be sufficient 
to determine their critical curves accurately. This is indeed 
the case in the largest k nown lens clust ers such as A1689 
(iBroadhurst et al.ll2005bh . CL002 4+1654 dKneib et al.l | 2003 ; 
Jee et al.ll2007h. RX J1347-1145 dBradac et al.ll2005l I2OO8 ; 
Halkola et al.ll2008l ). SDSS J1209+2640 <|Ofek et al.ll2008h . 
and RCS2 2327-02 (M. Gladders et al, in preparation) in 
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Figure 10. The largest Einstein radius in all-sky as a function 
of primordial non-Gaussianities /nl- WMAP5 cosmology is as- 
sumed. The error-bars show 68% confidence estimated from 300 
Monte-Carlo realisations. The source redshifts are assumed to be 
z s = 1, 3, and 7 from top to bottom panels. 



cal surveys, s uch as done by the Large Synoptic Survey Tele- 
scope (LSST; llvezic et al.lf2008l fl provide a promising way 
to locate clusters with very large Einstein radii. In such opti- 
cal surveys, we will be able to search for strongly lensed arcs 
directly without a priori information on locations of mas- 
sive c lusters (e.g.. iHoresh et al"1l2005l ; ISeidel fc Bartelmannl 
l2007f). Such blind/autom ated arc survey has been attempted 
by Cabanac et al.l (|2007T ) using Canada- France-Hawaii Tele- 
scope Legacy Survey (CFHTLS) data. The detection and 
characterisation of m assive clusters using weak lensing (e.g., 
iMivazaki et""aH I2007T ) will complement the identifications 
of lensed arcs and will allow a check of the model of 



which many strong lensing events are already identified. 
The discussion above suggests that wide-field deep opti- 
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Figure 11. Similar to Figure flOl but the numbers of clusters 
with the Einstein radii larger than certain values are plotted. 
The limiting Einstein radii are 20" (z s = 1), 35" (z s = 3), and 
40" (zs = 7). 



clusters that we adopted in Section [2] Another approach 
is to make use of (shallower) optical, X-ray, or Sunyaev- 
Zel'dovich (SZ) surveys to identify candidate massive clus- 
ters, and conduct follow-up optical imaging of each mas- 
sive cluster to characterise its lensing properties. Exam- 
ples of such optical/X-r ay/SZ cluster surve ys include the 
maxBCG cluster survey (|Koester et al.ll2007h . the ROSAT- 
ESO Flux Limited X- ray (REFLEX) Galaxy cluster survey 
llBohringer et alj|2004) . the Massive Cluster Survey (MACS; 
lEbeling et al.ll2007l) the RO SAT PSPC Galaxy Cluster Sur- 
vey dBurenin et al.ll2007t). th e Red-Sequence Cluster Survey 
(RCS; iGladders fc Yeell2005l l and planned SZ cluster sur- 
veys such as the South Pole Telescope (SPTfl the Atacama 



Cosmology Telescope (ACTn, the Atacama Pathfinder Ex- 
periment (APEX) SZ survejQ and a survey using the Planck 
satellite. However, clusters with largest Einstein radii are 
not necessarily the most massive. Cluster surveys need to 
be deep enough to locate masses as small as ~ 5 x 1O 14 M 
to assure completeness (see also Figure 2J. 

The critical curves of the largest lenses predicted in 
our model may offer guidance for identifying such systems 
in observations. In Figure 1121 we plot our prediction for 
the plausible critical curves of the larges t lens, as well as 
the crit ical curves of A1689 obtained in iBroadhurst et all 
l)2005br ). Because of the high concentration and rounder 
shape of the projected mass distribution, the predicted 
inner critical curves are rather small compared with the 
outer critical curve. Therefore for these systems strong lens 
events are dominated by standard "double" and "quadru- 
ple" image configurations. This is in marked contrast to 
typical clusters in which less concentrated 2D mass dis- 
tributions increase the importance of their inner criti- 
cal c urves and produce naked cu s p image configurations 
(e.g. iBlandford fc Kochanekl 1 19871 ; lOguri fc Keetonl |2004| ; 
iQguri et al.ll200ct ). We find that the critical curves of the 
largest known lens cluster A1689 are similar despite the per- 
turbation by several substructures. 

One of our main findings is that large lens clusters rep- 
resent a highly biased population. Whilst the biases in pro- 
jected 2D mass distributions can directly be tested from 
lensing observations of clusters, the alignment between ma- 
jor axes and line-of-sight directions will require additional 
observations. For instance, the line-of-sight elongation of a 
cluster can be inferred by combining multi-wavelength data 
such as X-ray, SZ , and kinematics o f member galax i es (e.g., 
Fox fc Pen! |2002l; iLee fc Sutol 12004 iGavazzil 120051 : ISerenol 



2007t lAmeglio et alj|2007[ ) 



6.2 Effect of baryons 

The triaxial halo model of JS02, which we adopted, 
is based on iV-body simulations of dark matter. It 
is of interest to see how baryon physics can af- 
fect our results. The most important baryonic ef- 
fect on clust er strong lensing comes from the c entra l 



feet on clust er str ong lensing comes from the c entra . 
galaxy (e.g.. iMeneghetti et all l2003bl ; ISand et all 120051: 



Puchwein et all 120051; iHennawi et al.l [2007bl: 



2008al ; Wambsganss et al.l l200St iHilbert et all 



Rozo et all 



200S 



Al- 
though the central galaxy can boost the lensing probabil- 
ity as high as ~ 100%, the effect is clearly scale dependent 
such that clusters with smaller Einstein radii are more sub- 
stantially affected by central galaxies. Both analytic and nu- 
merical studies agree in that for lensing with large Einstein 
radii (#e > 20") the enhancement of lensing probabilities 
due to cent r al galaxies become neglig i bly small (see, e.g. , 
lOguril |200d I Wamb sganss et all 12008 ; IHilbert et all I2OO8I ; 
iMinor fc Kaplinghatil2008l ). suggesting that the clusters dis- 
cussed in this paper, which have unusually large Einstein 
radii, are not affected by central galaxies. This is also ex- 
pected from the fact that their critical curves extend far be- 



http:/ /pole. uchicago.edu 



4 http:/ /www. physics. princeton.edu/act/ 

5 http:/ /bolo. berkeley.edu/apexsz/ 
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Figure 12. Critical curves of A1689 (z s = 3) constra ined from 
multiple strong lens systems ( Broad hurst et al.ll2005bl) are com- 
pared with the most plausible critical curves of the largest lens in 
our model, which is obtained from the median parameter values 
listed Tabled] 



yond central galaxies. Therefore we conclude that the effect 
of central galaxies is negligibly small for our results. 

Baryons also influence the shape of clusters. For in- 
stance, dissipative gas cooling re sults in more spherical dark 
haloes dKazantzidis et al1l2004h . In addition, the inclusion 
of hot gas com ponents slightly enhances the con centration 
of dark haloes (|Rasia et al.ll2004j : iLin et al.ll200fj ). although 
current studies are insufficient to quantify its statistical im- 
pact. 



less than #e considered in this paper, and the effect cancels 
out to first order. 

Some clusters have quite complicated morphology 
which cannot be described by a simple ellipsoid. An ex- 
treme example is the merger as seen in the bullet cluster 
l|Clowe et al.ll2006h . w hich turned out to have a large im- 
pact on arc statistics (jTorri et al.l I2004T ). One of the rea- 
sons for the large effect of mergers on arc cross sections is 
that complicated mass distributions tend to induce many 
cusps in caustics where prominent, long and thin arcs are 
preferentially formed. In contrast, the size of the Einstein 
radius is simply determined by the mass it encloses, and 
therefore the Einstein radius should be less sensitive to the 
complexity of the mass distribution. In addition, clusters 
which undergo merger events can be represented approxi- 
mately by triaxial haloes with very small axis ratios a/c. In 
this sense our calculation includes the effect of mergers. The 
more robust treatment of mergers will require calibration 
using high-resolution iV-body simulations. 

In our calculations we have ignored any chance projec- 
tion of multiple clusters along the line-of-sight. We make a 
rough estimate for the expected number of such events as 
follows. We consider following two situations: (1) superpo- 
sition of two massive haloes, and (2) a small halo on top 
of a massive halo. For case (1), we require that both haloes 
should have masses M > 7 x 1O 14 M0, because the superpo- 
sition of two clusters with these masses result in the total 
Einstein radius of the system large enough to affect our re- 
sults (#e > 30" for z s = 3). From the all-sky number of such 
clusters, iV c iu ~ 4, 400 for WMAP5, we compute the chance 
prob ability as N^ST 1 ^ ]! + uj{6)\ w 0.1 x [1 + u(6)\ < 1 
(e.g.. iBrodwin et al.ll2007l for the angular correlation func- 
tion ui{9)), much smaller than the predicted abundance 
of large lenses (e.g., N(> 35") w 35 for z s — 3). For 
case (2), we consider superpositions of massive haloes with 
M > 1 x 10 15 M Q and smaller haloes with M > 1 X 1O 14 M 
from the same reason, i.e., the superposition of these haloes 
with typical concentrations yields the total Einstein radius 
with its size comparable to that discussed in the paper. 
The all-sky numbers, 1,100 and 86,000 respectively, sug- 
gest the chance probability of ~ 0.5 x [1 + uj(8)] < 1, again 
smaller compared with the predicted abundance of super- 
lenses. Thus we expect that the effect of the chance projec- 
tion of multiple clusters is not significant. 

However, it has been argued that A1689 may 
have a possible seco nd ary peak along lin e-o f-sight 
llAndersson fc Madeiskil 12(0041 : lLokas et all 120061 : 
IKing fc Corles? 2007), which implies that multiple clusters 
might in practice have some impact on the statistics of 
superlenses. Ray-tracing in large-box TV-body simulations 
should provide an important cross-check of how important 
projections of haloes along line-of-sight are. 



6.3 Other possible systematics 

Cluster substructures affect the shapes and locations of 
arcs and hence m a y infl uence the Einstein radius. However, 
iMeneghetti et aD 1 2007h found t hat the radial shift due to 
substructures is < 5" (see also IPeirani et al.l 120081 ) , much 



7 CONCLUSION 

We have calculated the expected distributions of large Ein- 
stein radii in all-sky (40,000 deg 2 ) using a triaxial halo 
model. Our approach to generate all-sky mock catalogue of 
massive clusters and the properties of individual clusters 
with the Monte-Carlo method allows us to evaluate the cos- 
mic variance for such statistics, and at the same time, to 
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study biases in the population of clusters having such large 
Einstein radii. 

The largest Einstein radius in all-sky for source redshift 
z s = 1 was predicted to be 427_y arcseconds for WMAP5, 
35j_g arcseconds for WMAP3, and 54l^ 2 arcseconds for 
WMAP1, where errors are la cosmic variance. The sensi- 
tivity to as suggests that this statistic is a good measure of 
it. The Einstein radii are approximately twice as large for 
larger source redshift, z s — 7. In some realisations the largest 
lens cluster is the most massive cluster in the universe; in 
others smaller than 10 15 Mq. We have found that the popu- 
lation of these "superlens" clusters are significantly biased: 
their major axes are almost always aligned with the line- 
of-sight, and their projected 2D mass distributions appear 
rounder (by Ae ~ 0.2) and more concentrated (~ 40 — 60% 
larger values of concentration parameters). These biases are 
stronger than those fou nd in more common len s clusters with 
smaller Einstein radii (jHennawi et al.ll2007bh . In particular 
we have pointed out that the high concentration observed in 
A1689 is consistent with our theoretical expectation at the 
1.2a level. Thus the combined analysis of several clusters 
will be essential to address the claimed high concentration 
problem. Finally, we have studied the effect of primordial 
non-Gaussianities, and concluded that the abundance of rel- 
atively large lens clusters can in principle constrain primor- 
dial non-Gaussianities at a level comparable to the current 
CMB experiments (|/nl| ~ 100), if other cosmological pa- 
rameters are fixed. 

It will be very important to compare our analytic 
predictions with ray-tracing in TV-body simulations. In 
particular, the large cosmological Millennium Simulation 
l|Springel et al.l I2005T ) has sufficient resolution to resolve 
the centres of massive dark h aloes for strong lensing stud- 
ies (|Hilbert et alj|2007l . lioosh . Although its small box size 
(500/i _1 Mpc) does not allow predictions for all-sky distribu- 
tions of Einstein radii, we can use these simulations to val- 
idate and calibrate our semi-analytical model predictions. 
The comparison of our results with those from TV-body sim- 
ulations will be presented in a forthcoming paper. 

The statistics of large Einstein radii provide an impor- 
tant opportunity to test the standard FACDM paradigm, as 
it probes both the high-mass end of the cluster mass func- 
tion and central mass distributions of massive clusters. The 
measurement of Einstein radii is fairly robust, and future 
all-sky samples will soon be available to perform this study. 
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